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Abstract. The possibility of measuring the proton electromagnetic form factors in the time-like region at 
FAIR with the PANDA detector is discussed. Detailed simulations on signal efficiency for the annihilation 
of p + p into a lepton pair as well as for the most important background channels have been performed. It 
is shown that precise measurements of the differential cross section of the reaction p + p — > e~ + e + can 
be obtained in a wide kinematical range. The determination of the ratio 1Z of the moduli of the electric 
and magnetic proton form factors will be possible up to a value of momentum transfer squared of q 2 ~ 14 
(GeV/c) 2 with absolute precision from 0.01 to 0.5 (for 1Z ~ 1). The total p + p — > e~ + e + cross section 
will be measured up to q 2 ~ 28 (GeV/c) 2 . The results obtained from simulated events are compared to 
the existing data. Sensitivity to the two photons exchange mechanism is also investigated. 

PACS. 25.43.+t Antiproton-induced reactions - 13.40.Gp Electromagnetic form factors 



1 Introduction 

The availability of a high intensity antiproton beam up 
to a momentum of 15 GeV/c at the FAIR facility [I] and 
of the PANDA detector offers unique possibilities for new 
investigations in the field of hadron structure (see [2] for a 
review). Here we focus on feasibility studies for the de- 
termination of the proton electromagnetic form factors 
(FFs), in the time- like (TL) region [3], through the an- 
nihilation reaction: 



P + P 



(1) 



The underlying mechanism is assumed to be the exchange 
of one virtual photon of four momentum squared q 2 . The 
sensitivity of the measurement to higher exchanges, which 
are in principle suppressed, is also investigated. Muons 
carry the same physical information on the nucleon struc- 
ture as the electrons, but this work will focus on the elec- 
tron channel only. Although the measurements of electro- 
magnetic nucleon FFs have been going on since more than 
fifty years, major progress has been recently achieved in a 
wide kinematical region, mostly in space-like (SL) through 
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polarized elastic electron proton scattering. The determi- 
nation of form factors is limited by the steep decrease of 
the cross section with q 2 . Moreover in the time-like (TL) 
region measurements in both channels: p + p -f* (r + £ + 
are scarce and affected by poor statistics. 

The intensity of the antiproton beam, together with 
the performances of the PANDA detector, will make pos- 
sible the determination of FFs up to large q 2 . FFs are ex- 
tracted from the angular distribution of one of the charged 
lcptons. In reaction (JlJ, the difficulty of the measurement 
is related to the hadronic background, mostly annihilation 
into pions, which is six order of magnitudes larger than the 
production of a lepton pair. In this paper we report on de- 
tailed simulations of the hadronic background and discuss 
the precision and the significance of the extracted data in 
a wide kinematical range. 

This paper is structured as follows. In sect. [H the in- 
terest of measuring proton electromagnetic form factors 
is briefly recalled, and the present experimental situation 
both in SL and TL regions is illustrated. In sect. |3] simula- 
tion studies for both signal and background are reported. 
In sect.2]the results on the extracted FFs and on their pre- 
cision are discussed and compared with the existing data. 
The sensitivity to the two photon exchange mechanism 
is discussed in sect. [5] In sect. [5J the PANDA expected 
performances are compared to the ones obtained in previ- 
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ous experiments. The main results are summarized in the 
conclusions. 



2 Physics motivation 

Hadron electromagnetic FFs describe the internal struc- 
ture of a particle. Elastic FFs contain information on the 
hadron ground state, and are traditionally measured using 
electron hadron elastic scattering, assuming that the in- 
teraction occurs through one-photon exchange (OPE). As- 
suming a Parity and Time invariant theory, a hadron with 
spin 5 is described by 25 + 1 independent FFs. Protons 
and neutrons (spin 1/2 particles) are then characterized 
by two form factors, an electric Ge and a magnetic Gm 
which are analytical functions of one kinematical variable 



2.1 Space-like region 



the polarization method [9"lfTU] , show that the electric and 
magnetic distributions in the proton are different, contrary 
to what was earlier assumed. The Q 2 dependence of Ge 
and Gm, deduced from polarization experiments p(e,e)p 
differs from ([3]). The FF ratio shows a linear deviation 
from a constant, which can be parametrized as [TTj : 

(iGe/Gm = 1-059 - 0.143 Q 2 [(GeV/c) 2 ] for Q 2 > 0.4, 

(4) 

up to at least Q 2 =5.8 (GeV/c) 2 . Polarization measure- 
ments have recently been extended up to Q 2 =8.5 (GeV/c) 2 
by the GEP collaboration, at Jefferson Laboratory (JLab) 
PTj and may show a zero crossing for this ratio, if the lin- 
ear extrapolation of the fit, eq. (j4]) will be confirmed by 
the final results. 

As no experimental bias has been found in the ex- 
periments, the discrepancy between FFs determined from 
polarized and unpolarized measurements, has been at- 
tributed to radiative corrections, as two photon exchange 
(TPE) daHamHUITllin] or higher order corrections QI] . 



Elastic electron proton scattering allows to access the SL 2.2 Time-like region 



region, where FFs are real functions of Q 



-q 2 > 0. 



Electromagnetic FFs are determined through the e de- 
pendence of the (reduced) elastic differential cross section, 
which may be written, in OPE approximation, as [4]: 



,Q 2 ) = 



i 



2-|sm 2 (0 e /2) 



4E 2 sin> e /2) 



, da 

x < 1 + T) do 



a 2 cos 2 (6> e /2) 
rG 2 M (Q 2 )+eG 2 E (Q 2 ), (2) 



[1 + 2(1 + t) tan 2 (6> e /2)] 



Q 2 

UP ' 



a e = 1/137, is the electromagnetic fine structure constant, 
M is the proton mass, E is the incident electron energy 
and 9 e is the scattering angle of the outgoing electron. 
Measurements of a re d(0 e , Q 2 ) at different angles for a fixed 
value of Q 2 allow to extract Ge(Q 2 ) and Gm(Q 2 ) from 
the slope and the intercept of the linear e dependence (0) 
(Rosenbluth separation) . 

The existing data on G m are described by a dipole be- 
havior up to the highest measured value Q 2 ~ 31 (GeV/c) 2 
[5], according to: 



G M (Q 2 )/v = G d (Q 2 ), G d (Q 2 ) = (1 + Q70.71)" 



(3) 



where /j, ~ 2.79 is the proton anomalous magnetic moment 
in nuclear magnetons and Q 2 is expressed in (GeV/c) 2 . 
The independent determination of Gm and Ge from the 
unpolarized e~p cross section has been obtained up to 
Q 2 = 8.8 (GeV/c) 2 [5], and gives G E ~ Ga//m- Further 
extraction of Gm [5] is based on this assumption. 

Experimental and theoretical studies have been done 
since a few decades, but recent interest aroused due to the 
possibility to reach higher precision and larger values of 
q 2 at high intensity accelerators, using polarized beams, 
targets and polarimeters in the GeV range [TJ. In partic- 
ular recent measurements of the FF ratio [5], based on 



The TL region, where q 2 > 0, can be investigated using 
the crossed reactions p + p ■<-» e~ + e + . Due to unitar- 
ity, hadron FFs are complex functions of q 2 , and their full 
determination requires more observables as shown in |19l 
[20] . and recently discussed in 21 . However the unpolar- 
ized cross section depends only on their moduli, and their 
measurement is, in principle, simpler than in SL region. 
In SL region, the Rosenbluth separation requires at least 
two measurements at fixed q 2 and different angles, which 
implies a change of incident energy and scattered electron 
angle at each q 2 point. In TL region, the individual deter- 
mination of \Ge\ and \Gm\ requires the measurement of 
the angular distribution of the outgoing leptons, at fixed 
total energy s — q 2 . Previous experiments (see [22]), have 
measured the cross section up to q 2 =18 (GeV/c) 2 and ex- 
tracted \Gm\ in the hypothesis Ge = Gm or Ge = 
(which affects up to 30% the values of | G a/ | ) - 

Attempts to determine the ratio 1Z = |G_e|/|Ga:H can 
be found in the literature, in ref. [23] (PS 170 at LEAR) 
and more recently in ref. [24] . through measurements of 
the initial state radiation reaction (ISR) e~+e + — > p+p+"f 
(BABAR Collaboration). The results of the two exper- 
iments, although affected by large errors, seem to show 
a different trend. In the second case a larger value was 
found, in a wide q 2 range above threshold. 

The detector PANDA, using the antiproton beam plan- 
ned at FAIR, will open a new opportunity to measure TL 
FFs. The aim of this paper is to show the precision that 
can be achieved in the measurements of TL proton FFs at 
PANDA in a wide range of q 2 . The interest in spanning a 
large kinematical domain is the investigation of the transi- 
tion region from soft to hard scattering mechanisms, which 
is the domain of perturbative Quantum Chromodynamics 
(pQCD), where the nucleon can be described in terms of 
quark and gluon degrees of freedom. In such region, scal- 
ing laws and helicity conservation |25 1 I26 | give predictions 
for the asymptotic behavior of FFs. 
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Moreover, the comparison of SL and TL data, allows 
to verify asymptotic properties which hold for analytical 
functions [271128] . Following the Phragmen-Lindelof the- 
orem FFs in TL and SL region have to coincide 
for \q 2 \ — > oo. This implies not only that the moduli 
should be the same, a feature that will be tested with 
PANDA, but also that the phases of TL FFs, which can 
be accessed only through polarization measurements [TO1 
[201111], should vanish. 



where T> is obtained fitting the experimental data and 
A = 0.3 GeV is the QCD scale parameter. The calcu- 
lated cross section <j{<Jqcd) and the number of counts 
N(Nqcd) are given in table Q] assuming an integrated 
luminosity of C = 2 fb _1 , which is expected for each data 
point in four months data taking, with 100 % efficiency 
and full acceptanc^B It is assumed that \Ge\ = \Gm\, 
calculated from Eqs. flTJ and ©, respectively. The event 



3 Simulation studies 

3.1 Differential cross section and counting rate 

The differential cross section for the annihilation process 
([1]) , first obtained in ref. [T5] , is expressed as a function of 
the proton electromagnetic FFs as: 



do 



d(cos9) 



8M 2 t^/t(t - 1) 

'E 



[r|G M | 2 (l + cos 2 



|G B | 2 sin 2 l 



(5) 



where is the electron production angle in the center 
of mass system (CM). The cos 2 9 dependence of eq. ([5]) 
results directly from the assumption of OPE, where the 
spin of the photon is equal to one and the electromag- 
netic hadron interaction satisfies C invariance. This cor- 
responds, by crossing symmetry, to the linear Rosenbluth 
cot 2 (6> e /2) dependence [30]. 
The total cross section is: 



7TOC 



6M 2 Ty/r(T - 1) 



2t\G 



M 



(0) 



The evaluations of the cross section and of the counting 
rate require the knowledge of the FFs. For the numerical 
estimates below, we use a parameterization of |Gj\/| from 
l27l , where the numerator is a constant fitted on TL data: 



\G M \ = 22.5 (1 + q 2 /0.7l) 2 (1 + g 2 /3-6)" 



(7) 



Here q 2 is expressed in (GeV/c) 2 . Eq. (0 gives a con- 
servative estimation of the yield at large q 2 . As the TL 
\Gm\ values have been extracted from cross section mea- 
surements assuming \Ge \ = \Gm\, the same hypothesis is 
taken for counting rate estimates, on the basis of eq. (jOJ. 
The evaluations of the cross section and of the counting 
rate have been also performed using the following QCD 
inspired parameterization of \Ge,m\, based on analytical 
extension of the dipole formula eq. ([3]) in TL region, where 
Q 2 is replaced by q 2 . Corrections based on dispersion re- 
lations have been suggested in [3T] to avoid 'ghost' poles 
in a s (the strong interaction running constant), and can 
be included in the following form: 



Iff 



QCD, 
E,M I 



V 



[log 2 {s/A 2 ) 



V = 89.45 [GeV/c] 



s 


V 




N 


OQCD 


Nqcd 


[GeV/c] 2 


[GeV/c] 


[Pb] 




[Pb] 




5.40 


1.7 


538 


1.1 10 6 


481 


9.6 10 5 


7.27 


2.78 


72 


1.4 10 5 


69 


1.4 10 5 


8.21 


3.3 


32 


6.4 10 4 


33 


6.5 10 4 


11.0 


4.9 


4.52 


9.1 10 3 


5.48 


1.1 10 4 


12.9 


5.9 


1.6 


3.2 10 3 


2 


4.3 10 3 


13.8 


6.4 


1 


2 10 3 


1.4 


2.8 10 3 


16.7 


7.9 


0.29 


580 


0.49 


979 


22.3 


10.9 


0.04 


81 


0.09 


183 


27.9 


13.4 


0.01 


18 


0.03 


51 



Table 1. Cross section a (ctqcd) and number of counts, N 
(Nqcd) from eq. [7] (eq. [8]) corresponding to an integrated lu- 
minosity of £ = 2 fb _1 , for different values of q 2 = s and of 
the antiproton momentum, p. 



generator for the reaction (fT]), is based on the angular 
distributions from eq. (|5"]). with prescription ([Jj for the 
magnetic form factor Gm- 

Three different hypothesis were taken for Ge- Besides 
the case \Ge\ = \Gm\, (TZ = 1), which is strictly valid 
only at threshold, the case TZ = and the case TZ = 3 
(as suggested in ref. [32]), were also considered. The cor- 
responding angular distributions were built keeping the 
same total cross section at each q 2 . They are shown in fig. 
CD for three values of q 2 = 5.4, 8.2, 13.8 (GeV/c) 2 . The 
reported error bars are statistical only. The sensitivity to 
TZ decreases when q 2 increases, due the falling of the cross 
section and to the relative weight of the magnetic term, 
which is growing as q 2 . 



3.2 Detector description 

An extensive description of the PANDA detector and of 
its different components can be found in ref. [2]. In the 
present work, we mention some of the characteristics which 
are important for the specific processes of interest. The 
conception of the detector, the read out and the acquisi- 
tion benefit from the experience gained from the construc- 
tion of recent detectors, such as ATLAS, CMS, COMPASS 
and BABAR. The high quality antiproton beam of mo- 
mentum from 1.5 to 15 GeV/c will be provided by the 
High Energy Storage Ring (HESR), equipped with elec- 
tron and stochastic cooling systems. The PANDA detector 



(8) 



1 This value will be always used below, except when ex- 
plicitely indicated. 
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Fig. 1. (Color online) CM angular distributions from the event 
generator, at q 2 —5A, 8.2, and 13.8 (GeV/c) 2 , for p+p —¥ e~e + 
and three different hypothesis: 1Z = (black solid circles), 
TZ = 1 (red triangles), and TZ = 3 (black open circles), keeping 
the same value of the total cross section. 



should ensure the detection of 2-10 7 interactions/s, with 
47T acceptance and momentum resolution for charged par- 
ticles at a few percent level. The expected average lumi- 
nosity L = 1.6 • 10 32 cm -2 s _1 will be reached with a 
pellet target of thickness 4 • 10 15 hydrogen atoms/cm 2 , 
and 10 11 stored antiprotons in HESR. The target will 
be surrounded by a spectrometer which includes a 2 T 
solenoid magnet. A forward spectrometer, based on a 2 
Tm dipole magnet, ensures the detection at small angles, 
down to 2°. Each spectrometer is equipped with detectors 
for tracking, charged particle identification, electromag- 
netic calorimetry and muon identification. The interac- 
tion point is reconstructed with the help of a microver- 
tex detector (MVD) , consisting of layers of radiation hard 
silicon pixel detectors surrounded by silicon strip detec- 
tors. In addition, charged particle tracking and identifi- 
cation is provided by straw tubes (STT) or a time pro- 
jection chamber (TPC), completed by GEM detectors at 
forward angles. The identification of hadrons and leptons 
in a wide kinematical range requires other complementary 
detectors. The time-of-flight of particles emitted at large 
polar angles will be measured in a good resolution time- 
of-flight barrel. The electromagnetic calorimeter (based on 
PbW04 crystals) will provide good energy and time reso- 
lution for the detection of photons and electrons at inter- 
mediate energy from a few MeV to ~ 10 GeV. Detectors 
based on Cherenkov light (DIRC), which are very efficient 
for pion-electron separation for momentum p < 1 GeV/c, 
will be used in a barrel detector and a forward endcap 
detector. 



3.3 Simulation, digitization and reconstruction 

The simulations are based on the same PANDA software 
and the same detector geometry as the one used for the 
physics benchmark simulations presented in [5] , with STT 
as central tracker and a pellet target. The simulation con- 
sists of two steps. The first one, based on the GEANT4 
code, is the propagation of the particles through the de- 
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Fig. 2. (Color online) Ratio of the energy deposit in the elec- 
tromagnetic calorimeter to the particle momentum E/p as a 
function of p, for pions and electrons separation. 



tector. The information on the hits and the energy losses 
has been digitized, including a model for electronic noise, 
into a response of the different detectors, in view of using 
the same Monte Carlo code in simulations and in future 
real data analysis. 

The second step is the reconstruction of the physi- 
cal quantities particularly important for electron identi- 
fication such as momentum, ratio of energy loss to path 
length in each straw tube (dE/dx), Cerenkov angle in the 
DIRC detector, and energy deposit in the electromagnetic 
calorimeter. 

These two steps have been described in detail in [2] 
and we will give here only the details which are the most 
important for the electron identification. 

The parameters used to simulate the fluctuations of 
the physical signals have a direct influence on the parti- 
cle identification capabilities. Precise estimate of energy 
loss fluctuations in thin layers are needed for the central 
tracker, and were included according to the PAI model for 
the description of the ionization process |33j . 

The dE/dx values are then used for particle identifica- 
tion using the truncated arithmetic mean method in order 
to exclude from the sample the largest values correspond- 
ing to the extended Landau tail. A truncation parame- 
ter corresponding to 70% out of the N individual dE/dx 
values was taken to calculate the arithmetic mean, as a 
compromise between the requirements of the best reso- 
lution, defined as the width of the gaussian fit, and the 
smallest tail of the distribution. A resolution of < 10% 
is obtained for pions at 1 GeV/c, which corresponds in 
average to a value of 4cr of the distance between the two 
truncated means for electrons and pions. For the DIRC, 
the resolution on the Cerenkov angle is given by 
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with the single photon resolution oc^ = 10 mrad. The 
number of detected photons N p h depends on the velocity 
and path length of the particle within the radiator and 
takes into account transmission and reflectivity losses as 
well as the quantum efficiency of the photodetectors. A 
nearly Gaussian resolution of about 2.3 mrad is obtained 
for 1 GeV/c pions [34]. As the Cerenkov angles for pions 
and electrons differ by 36 mrad at 500 MeV/c and by 4 
mrad at 1.5 GeV/c, the DIRC has a significant discrimi- 
nation power at the lowest energies. 

The electromagnetic calorimeter is the most important 
detector for the electron identification through the ratio 
E/p of the measured energy deposit to the reconstructed 
momentum (see fig. [2]). Electrons deposit all the energy 
in an electromagnetic shower, while muons and hadrons 
loose only a much lower fraction of their kinetic energy by 
ionization processes. However, high energy deposits may 
be due to hadronic interactions within the crystals. In par- 
ticular, for charged pions undergoing quasi elastic charge 
exchange processes almost all the energy is transferred to 
a neutral pion decaying into two photons. These processes 
are taken into account by choosing a GEANT4 physics list 
which includes the Bertini intra-nuclear Cascade model for 
hadron interactions at intermediate nuclear energies from 
hundreds of MeV to GeV [35]. 

The shower shape can also be used for particle identifi- 
cation. Due to the small Moliere radius (2 cm) of PbW04, 
of the order of the crystal size (2.2 cm), the largest fraction 
of the electromagnetic shower is contained in a few mod- 
ules, while a hadronic shower with similar energy is more 
spread. The shower shape is characterized by the energy 
deposits in the central crystal and in the 3x3 and 5x5 mod- 
ule arrays containing the central scintillator. In addition, 
a set of four Zernike moments [3 6) describes the energy 
distribution within the shower by polynomials which are 
functions of radial and angular coordinates. 

Sets of particles of different species have been gener- 
ated at given momenta and polar angle. The probabili- 
ties for identification of a given particle as electron, muon 
pion, kaon, or proton are then calculated for each detec- 
tor using, in addition to the variables discussed above, 
the dE/dx information in the microvertex detector, and 
the information from the muon detector. For the electro- 
magnetic calorimeter, the probability is calculated from 
the output of a neural network using as inputs all the pa- 
rameters of the shower listed above, as described in [2J. 
As an example, the probability for a pion to be identi- 
fied as an electron is of the order of 10 -3 at p=3 GeV/c, 
which agrees with values measured on existing detectors 
(BABAR). From the individual subdetector likelihoods, 
a global PID likelihood is then calculated. Depending on 
the signal and background channels, the thresholds can 
be adjusted in order to ensure the required purity while 
keeping the signal efficiency as high as possible. In our 
case, as the pion yield is much larger than the electron 
one, the threshold is defined by the purity requirement on 
the electron signal. 



3.4 Study of background channels 

Reactions involving two or more hadrons in the final chan- 
nel constitute an important background for the measure- 
ment of channels with creation of a lepton pair. 

Due to the hermiticity of the detector and the good 
tracking resolution, channels with three or more hadrons 
in the final state will be very efficiently identified. The 
cross section for channels involving three pions is known 
to be at most an order of magnitude larger than two pion 
production. Our simulation show that a reasonable cut on 
the missing mass gives a rejection factor of at least a factor 
of hundred. 

The larger background is then expected to come from 
annihilation into two hadrons &s p + p — > tt° + tt°, p + p — > 
7r~ +7T+ oi p + p — > K~ + K + . The cross sections for the 
neutral (charged) channels production are about five (six) 
orders of magnitude larger than for reaction ([T]). In the 
case of the ir° + tt° production, e~e + pairs are produced 
after conversion of the photons from the main tt° decay, 
in particular in the beam pipe before the tracking system. 
In addition, one (or both) ir° may undergo Dalitz decay, 
7T° -)• e~ + e+ + 7, with probability 10~ 2 (10~ 4 ). 

In case of charged hadron pair production, both hadrons 
can be misidentified as leptons. In case of kaon production, 
the probability of misidentification is lower and kinemati- 
cal constraints are more efficient, due to their larger mass. 

Therefore, the background coming from pp annihila- 
tion into two pions is expected to be the largest and has 
been evaluated using detailed simulations. 

3.4.1 Simulations of p + p — > tt~ + ir + and 
p + p — > 7T° + 7T° reactions 

The angular distributions for charged and neutral pion 
pair production were extrapolated from a parameteriza- 
tion of the data L37 , 38 , 39 , ^SUS2] - For s < 6 (GeV/c) 2 , 
the existing data [37ll40ll4T] were fitted by Legendre poly- 
nomials. In the high energy range, instead, the behavior 
of exclusive processes is driven by dimensional counting 
rules, thus the differential cross section of the p + p — >• 
7r~ + 7T + process can be parametrized as [231126] : 



where 9 is the CM angle of the 7r~, t is the Mandelstam 
variable and the function f(9) depends on the reaction 
mechanism. In the framework of the quark interchange 
dominance model |43j . one has 

f(0) = ±(l-z 2 )[2(l-z)- 2 + (l + z)- 2 } 2 , z = cos9. (10) 

C = 440 mb (GeV/c) 14 is a constant, which can not 
be predicted by QCD, and it is determined from ir + p 
elastic scattering at momentum 10 GeV/c and cos 9=0. 
The model predictions were symmetrized (da(9)/dt — > 
[da(9)/dt + d<j(n — 9)/dt}/2) and readjusted in the region 
around 90°, at each s value, to get a better agreement with 
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Fig. 3. (Color online) CM angular distributions used to sim- 
ulate the p + p — > 7r~ + 7r + reaction for s—5A (GeV/c) 2 (blue 
dotted line) for s=8.21 (GeV/c) 2 (green solid line) and s=13.5 
(GeV/c) 2 (red dashed line) as a function of the cosine of the ir 
production angle. Data are from ref. ,37. (squares) and from 
ref. [381139] (circles). 



the data. The results of the event generator are shown in 
fig. ED for s=5.4 (GeV/c) 2 , s=8.21 (GeV/c) 2 and s=13.5 
(GeV/c) 2 and compared to data obtained in refs. [27,38, 

129]. 

For exclusive 7r°7r° production at high energy the fol- 
lowing parameterization was taken: 

^ = J^W /(•••J-XX-we-o ("» 

where Pi(cos0) are Legendre polynomials and are fitted 
to the data from E760 at Fermilab, in the kinematical 
range 8.5 < s < 18.3 GeV/c as shown on fig. gj The 
quality of the fits can be seen in fig. (Q| , where examples 
of differential cross sections for p+p -+ ir° + 7r° are shown. 

For both reactions, the extrapolation to | cos6>| = I is 
affected by a large uncertainty, as no data exist at high en- 
ergy. Therefore, in the following, only the angular region 
| cos#| < 0.8 will be considered, in particular for the eval- 
uation of statistical and systematic error sfl In this region, 
the ratio of 7r~ + ir + to e~ + e + cross sections varies from 
I0 5 at | cos6>| = to 3-10 6 at | cos6»| = 0.8. On the basis of 
these numbers, the rejection power should be larger than 
3 • 10 9 (3 • 10 8 ) for 7r~ + 7r + (ir° + 7T°) in this angular range 
to limit the background to 0.1% level. 

In order to evaluate the background rate fulfilling the 
e~ + e + criteria, the generated ir~ +7r H 
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Fig. 4. (Color online) Angular distribution of pions produced 
in p + p — >■ 7T° + vr° annihilation, for s = 8.47 (GeV/c) 2 (red 
circles) and s — 18.26 (GeV/c) 2 (blue squares). Data are from 
ref. \M- 



were analyzed using the same PID cuts and kinematical 
constraints as for the signal. 



3.4.2 Suppression of ir ir + and 7r°7r° backgrounds 

Due to the difference of six order of magnitude in the cross 
section, between the signal and the 7r~7r + background, an 
event sample of at least 10 8 7r~7r + events was simulated 
at (? 2 =8.21, 13.8, and 16.7 (GeV/c) 2 . 

To discriminate pions from electrons, cuts have been 
applied to the PID combined likelihood for the assump- 
tion that the detected particle is an electron. The num- 
bers of simulated 7r _ 7r + events left after Loose, Tight and 
Very Tight PID cuts (corresponding respectively to mini- 
mum values of the identification probability 85%, 99% and 
99.8% for each lepton of the pair) are displayed in tabled] 
It is shown that the Very Tight cuts are needed to reach 
a rejection factor of a few 10 7 . 

Further selection based on the reaction vertex and on 
the kinematical fit method is applied. The kinematical 
fit method is a constrained fit, which takes into account 
energy and momentum conservation. From this, a confi- 
dence level (CL) associated to 7r + 7r~ hypothesis, CL^ and 
a second one, CL e corresponding to e + e~ hypothesis are 
calculated. The selection of the electrons results from two 
conditions: CL T . e > 10~ 3 (which corresponds to x 2 < 7 
for the kinematical fit) and CL e > 10 CL^ . These condi- 
tions result in an additional rejection factor of the back- 
ground of ~ 100. Finally, combining the kinematical fit 
with the PID ends up in a overall background suppression 
factor of the order of a few 10 9 . 

It has been checked that, for | cos(0)| < 0.8, the ir~ir + 
contamination does not depend drastically on angle and 
will remain below 0.1% in the q 2 range of interest. Concerning 



the 7r°7r° channel, Dalitz decay, 7r° 



+ e+ + 7, has a 
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<? [GeV/c]^ 


8.2 


12.9 


16.7 


no cut 


10 s 


10 s 


2-10 s 




Loose 


425 


1.2-10 3 


3-10 3 


PID cuts 


Tight 


31 


70 


120 




Very Tight 


2 


5 


6 


kinematic fit(CL) 


8-10 5 


10 6 


2.5-10 6 



Table 2. Number of ir~iv + events, misidentified as e~e + , left 
after Loose, Tight and Very Tight PID cuts corresponding to 
respective minimum values of the electron identification proba- 
bility 85%, 99% and 99.8% and after the confidence level (CL) 
cut on the kinematic fit for three different q 2 values (see text). 
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Fig. 5. (Color online) Angular distribution of electrons from 
e~e pairs as a function of cos 6 at q 2 = 8.21 (GeV/c) 2 : gener- 
ated events (black circles), reconstructed events (red squares), 
acceptance and efficiency correction (blue stars, left scale), and 
efficiency corrected events (green triangles). 



be measured with PANDA is given in fig. [SJ for q 2 = 8.21 
(GcV/c) 2 , and assuming \Ge\ = \Gm\- The reconstructed 
events (red squares) are obtained after full Monte Carlo 
simulation which takes into account tracking, detector effi- 
ciency, and acceptance as described in sec. 13.31 The recon- 
struction efficiency corrections have been obtained from an 
independent simulation which assumes an isotropic elec- 
tron distribution (blue stars, right scale). Once corrected 
for this efficiency, the distribution (green triangles) nicely 
agrees with the generated one (black circles) . One can see 
that at this q 2 value the average efficiency is of the order 
of 40%. 

The reconstruction efficiency depends on the angle. 
The sudden drop at |cost9| > 0.8 corresponds to a de- 
crease of the PID efficiency. At g 2 =8.21 (GeV/c) 2 the poor 
dE/dx identification from the STT is responsible of this 
drop. The loss of efficiency at cos 6 = is due to the target 
system. 

The reconstruction efficiency, after integration over the 
angular range | cos#| < 0.8, is shown in fig. [5] It is maxi- 
mum at q 2 ~ 8 (GeV/c) 2 and decreases to 15% at q 2 ~ 23 
(GeV/c) 2 . The effects of PID and kinematical constraints 
are shown separately. The drop at large q 2 is mainly due 
to PID cuts, as the laboratory angular distribution is more 
forward peaked with increasing q 2 , whereas the kinemat- 
ical selection shows a rather constant behavior. 

The normalization of the measured counting rates will 
be provided using the pp luminosity detector, with an ex- 
pected precision of 3%. 

Standard radiative corrections are included in the sim- 
ulation program, via the PHOTOS package [H]. The data 
will have to be corrected for soft and hard photon emis- 
sion, which partially compensate each other, reducing the 
overall effect. Radiative corrections strongly depend on 
the kinematical conditions and on the criteria for data 
selection. At q 2 =9 GeV 2 rough estimate gives an overall 
effect of 10-15%. Detailed studies of the effects of radia- 
tive corrections for this specific channel will be object of 
a separate paper. 



probability 10~ 2 . Three processes can be sources of e~e + 
pairs: i) double Dalitz decay of the two w°, ii) Dalitz decay 
of one of the pion associated with gamma conversion from 
the other pion, iii) photon conversion from two different 
pions. All these processes, with comparable rates, produce 
a six particle final state. Thus, even if the produced e~e + j 
pairs fulfill the PID cut, the kinematical constraints give 
a rejection factor, which combined with the 10~ 4 proba- 
bility for such processes, lead to an efficient suppression of 
this background. Moreover, by requiring that only a single 
e~ +e + pair has been identified in the whole detector solid 
angle, it is possible to reduce even further the contribution 
of this channel. 

3.5 Analysis of the e + e + channel 

Simulations were done for reaction ([1]) at q 2 values corre- 
sponding to tablcQ] for TZ =0, 1, and 3. A realistic projec- 
tion of the angular distribution of e~e + events, as it will 



4 Results and discussion 

For each q 2 value, the simulated differential cross section 
was fitted with a two-parameter function, in order to ex- 
tract a global normalization a and the form factor ratio 
TZ, according to: 

N(cos6) = a[r(l + cos 2 (9) + TZ 2 sin 2 9}. (12) 

The results are shown in fig. where the expected statis- 
tical uncertainty on TZ is plotted as a function of q 2 as a 
yellow band for the case 1Z — 1 , and compared with the ex- 
isting values from refs. [33] (squares) and [51] (triangles). 
Different methods were used to extract the errors bars for 
1Z. Different kinds of fits where performed with MINUIT, 
extracting TZ or TZ 2 from the quadratic expression (|12[) and 
also from the angular asymmetry which enters linearly in a 
cos 2 9 distribution (see sect. 5). A method based on Mon- 
tecarlo was also used. Detailed discussion and numerical 
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Fig. 6. (Color online) Reaction p + p —¥ e + + e~: average 
reconstruction efficiency as a function of q 2 (black solid line). 
The effect of applying kinematical constraints (red dashed line) 
and PID cuts (blue dotted line) are separately shown. 
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Fig. 7. (Color online) Expected statistical precision on the 
determination of the ratio 1Z, (yellow band) for 1Z = 1, as a 
function of q 2 , compared with the existing data from Refs. [23] 
(triangles) and [23] (squares). Curves are theoretical predic- 
tions (see text). 



values for the different values of 7Z can be found in ref. 
[4"5] . As an example, for g 2 =13.84 (GeV/c) 2 , one obtains 
from a Montecarlo approach: 1Z = 11q!51j ^ = ^iol53 an d 
an upper limit of 0.61 (CL=68%) for 1Z = 0. Therefore we 
concluded that a meaningful value for 1Z can be extracted 
up to q 2 ~ 14 (GeV/c) 2 . In the low g 2 region, the precision 
is at least an order of magnitude better than for the exist- 
ing data. With a precise measurement of the luminosity, 
this will allow to extract the moduli of Ge and Gm, to 
be compared with the corresponding SL values and with 
model calculations. 

Model predictions display a quite large dispersion, as 
shown in fig. A QCD inspired parameterization, based 
on scaling laws [251[26] , predicts 1Z = 1 , as it depends only 
on the number of constituent quarks (red dashed line). 
The green solid line is based on the vector meson dom- 
inance (VDM) approach from ref. [¥Bj, and grows up to 
q 2 ~ 15 (GeV/c) 2 . The blue dash-dotted line is the pre- 
diction from ref. [17], based also on VDM, but including 
terms to ensure the proper asymptotic behavior predicted 
by QCD. These models, originally built in the SL region, 
have been analytically extended to the TL region and the 
parameters have been readjusted in ref. |21j in order to 
fit the world data in the whole kinematical region (i.e., 
in SL region, the electric and magnetic proton and neu- 
tron FFs, and in TL region, the magnetic FF of the proton 
and the few existing data for neutron [IE])- Although these 
models reproduce reasonably well the FFs data, they give 
very different predictions for the form factor ratio. It is 
also shown in ref. [21] that polarization observables show 
large sensitivity to these models. 



At larger q 2 , when the sensitivity of the experiment 
will make the extraction of 1Z meaningless, it will then be 
possible to extract \Gm \ under a definite hypothesis on the 
ratio, in general TZ—1, as done in previous measurements. 
With a precise knowledge of the luminosity, the absolute 
cross section can be measured up to q 2 ~ 28 (GeV/c) 2 . 
The precision of such measurement is shown in fig. [5] The 
comparison with the world data shows an expected im- 
provement of at least a factor of ten. Here only the statis- 
tical accuracy, based on the number of events measured 
and identified is taken into account. The reported error 
bars are based on conservative extrapolation of the values 
reported in table 

Systematic effects of the tracking and reconstruction 
procedure will be evaluated on the real data, mostly by 
measurements on known reactions, which are the source 
of well controlled data samples. 

Comparing such data samples to simulations, as a func- 
tion of momentum and angle, will allow to check the rejec- 
tion power against pions for each subdetector and to de- 
termine the electron identification efficiency. This insures 
the control of the global rejection (purity of the electron 
data sample) and the knowledge of the overall electron 
reconstruction efficiency. 

The electromagnetic calorimeter is a crucial ingredi- 
ent for the electron/pion separation. Besides the calibra- 
tion using cosmic rays, several reactions can be identi- 
fied. The reaction p + p — > ir° + tt° can be used for the 
calibration of all electromagnetic calorimeter crystals at 
the sub percent level. The two body hadronic reaction 
p+p — > 7r + +7r~ allows to tune the PID parameters related 
to pions, and to check the calibration and the resolution 
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Fig. 8. (Color online) q 2 dependence of the world data on the 
effective proton TL FF, \Gm\, as extracted from the annihi- 
lation cross section assuming \Ge\ = |Gm|: BABAR [24] (red 
full circles); Fenice 48 (blue open circles ); E835 [22II49] (green 
open lozenge); PS 170 [23] (gray open stars ); E760 [42] (blue 
asterisk ); DM1 [50] (green full triangles); DM2 [ST] (green 
open squares); BES [52] (cyan open cross ); CLEO 53 (blue 
triangle down ) ; PANDA (full black squares) corresponding to 
an integrated luminosity of 2 fb , errors are statistical only 
(this work). 



of the momentum reconstruction at permille level, due to 
the high statistics. Systematic effects from detector mis- 
alignments are expected to be negligible. In addition, the 
decay J/W — > e + + e~ , with 6% branching ratio, provides 
electron data samples of good purity and known angular 
distribution. It will constitute a concrete measurement of 
the individual detector responses to electrons. At q 2 —9.6 
(GeV/c) 2 , one week measurement will be affected by 1% 
statistical precision in a 0.1 cos 8 bin. 



5 Sensitivity to two photon exchange 

As stressed in the introduction, the expression of the cross 
section (0) assumes OPE. TPE is suppressed by a factor 
of a e . At large q 2 , however, TPE could play a role, due in 
particular to a possible enhancement from a mechanism 
where the momentum is equally shared between the two 
photons 54,55,56,57 . Recently, the possibility of a siz- 
able TPE contribution has been discussed in connection 
with discrepancies between experimental data, on elastic 
electron deuteron scattering [30] taken in different exper- 
iments, and elastic electron proton scattering in polarized 
and unpolarized experiments [71[ni[T^[ni[TCl[Tl][T71[51[5g 
60 , 61j . Experimentally, no model independent and unam- 
biguous evidence of TPE (more exactly, of the real part of 
the interference between OPE and TPE) has been found 



in the experimental data ('>'_' li.'t (i 1 65^ . Recent analysis of 
e^p cross sections are controversial due to the uncertain- 
ties of the data |661|67J|68] . 

The general analysis of experimental observables in the 
reaction p+p — > e~+e + [69] and in the time reversed chan- 
nel [70] . taking into account the TPE contribution, was 
derived in a model independent formalism developed for 
elastic electron proton scattering [58,59,60 . It was shown 
that in presence of TPE, the matrix element contains three 
complex amplitudes: Ge, Gm and which depend on 
two kinematical variables, and that the angular distribu- 
tion contains new terms which are odd in cos 8 and are 
of the order of a e compared to the dominant contribution 

Therefore, it seems interesting to study the possibility 
of identifying the TPE effect in the present experiment. 
The purpose of this study is not to determine the physical 
amplitudes, but to set a limit for a detectable odd cos 9 
contribution, eventually present in the data. As the TPE 
amplitudes are not known, we used the presence of odd 
terms in cos 8 as a (model independent) signature, intro- 
ducing drastic approximations: we neglected those contri- 
butions to Ge,m which are smaller by an order a e . We 
considered only the real part of the three amplitudes, de- 
noted Ge, Gm and F 3 , as their relative phases are not 
known. 

We approximated the differential cross section in the 
following way: 



by taking 



dQ 4<? 2 V t - 1 



G 



(13) 



D ~ G' 2 M (1 + cos 2 6>) + — sin 2 8 + 

T 

2y/r(T- 1) [ — - Gm) F 3 cos 8 sin 2 I 



(14) 



Angular distributions were simulated according to cq. (1141) , 
assuming G E = G M for q 2 =5.4, 8.2 and 13.8 (GeV/c) 2 
and for F 3 /G M =2%, 5% and 20%. A number of events 
corresponding to Table I were processed for each of these 
q 2 values and the detector efficiency was taken into ac- 
count. The TPE components induce a distorsion in the 
angular distributions which vanishes at cos 8=0 and ± 1. 
In order to analyze the distributions, and extract the val- 
ues of the two photon amplitude, we rewrite the angular 
distribution as a polynomial in cos 8. In case of OPE, Eq. 
([5]) can be rewritten as 



da 



d(cos8) 



cr [l + A cos 2 



(15) 



where Co is the value of the differential cross section at 
8 = tt/2 and A is the angular asymmetry [27) . 

Therefore, at each value of q 2 , we can fit the angular- 
distributions with a straight line in cos 2 8, 



y = ao + ai cos 8, 



(16) 
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Fig. 9. (Color online) Angular distribution as function of 
cos 2 9, according to eq. HI]), for q 2 =5.4 (GeV/c) 2 and for 
different contribution of TPE : no contribution (top left), 2% 
contribution (top right), 5% contribution (bottom left), and 
20% contribution (bottom right). The solid (dashed) line is 
the result of the fit from eq. (|17fl (from eq. (116I) ) . Errors are 
statistical only. 



necessary to bring the data in agreement in the Q 2 range 
between 1 and 6 (GeV/c) 2 7 . The PANDA simulations 
show that such level of contribution will be detectable in 
the annihilation data. We have shown the stability of the 
extraction of ao and 02, from the data, even in presence of 
a relatively large contribution of TPE, in the approxima- 
tion (|14l) . This is due to the symmetry properties of the an- 
gular distribution and it has to be taken with caution: the 
relations between the observables (differential cross sec- 
tion and angular asymmetry) and FFs hold only in frame 
of OPE. The signification of the extracted parameters in 
terms of the moduli of the two electromagnetic FFs is not 
valid anymore. 

Let us stress that the main advantage of the search 
of TPE in TL region is that the information is fully con- 
tained in the angular distribution (which is equivalent to 
the charge asymmetry). In the same measurement, the odd 
terms corresponding to TPE can be singled out, whereas 
in SL region, in case of TPE, it is necessary to measure 
electron and positron scattering, in the same kinematical 
conditions. TPE effects cancel if one does not measure the 
charge of the outgoing lepton, or in the sum of the cross 
section at complementary angles, allowing to extract the 
moduli of the true FFs [591170]. 



6 Comparison with previous experiments 



where ao and a\ are related to the physical FFs. Devia- 
tions from a straight line are the evidence of the presence 
of higher order terms, beyond Born approximation. In or- 
der to check the sensitivity to odd terms, we fit the angular 
distributions by the function: 

y = a n + a% cos 2 9 + a 2 cos 0(1 — cos 2 9), (17) 

where a 2 is directly related to the ratio F^/Gm- 

The results of the fit are reported in table[3J In the case 
of OPE, as expected, the coefficient a 2 is compatible with 
zero. The odd cos 9 contribution starts to be visible for 
Fz/Gm > 5%. Note however that the extraction of 1Z and 
A is not affected, in the limit of the error bars, by the pres- 
ence of the C-odd term and that the ao and ai terms are 
very stable, even at large q 2 , although the statistical er- 
rors are more sizable. Here the error on 1Z is derived from 
A by first order derivation. Figure [9] shows the angular 
distribution as a function of cos 2 9, for g 2 =5.4 (GeV/c) 2 . 
The lower - red dots (upper - green dots) branches cor- 
respond to backward (forward) emission for a negative 
lepton. The solid (dashed) line is the result of the fit from 
eq. (fTT|) (from eq. ([To]) ), which includes (does not include) 
the odd cos 9 terms. 

Due to crossing symmetry properties, the reaction mech- 
anism should be the same in SL and TL regions, at similar 
values of the transferred momentum. If TPE is the reason 
of the discrepancy between the polarized and unpolarized 
FFs measurements in SL region, a contribution of 5% is 



As it has been stressed above, the HESR ring will pro- 
vide a high intensity antiproton beam. This feature, to- 
gether with the high performance of the PANDA detec- 
tor, will allow to have the best measurement until now of 
FFs in TL region. In table 2] we summarize kinematical 
and technical aspects of the existing TL FFs experiments. 
All previous results have been limited by low statistics, 
which prevented a precise determination of angular distri- 
butions. 

In the case of PS170 [33] the detector acceptance was 
limited, in particular due to the covering in azimuthal an- 
gle. The large efficiency for E835 [49] is due to the detec- 
tion based on gas Cerenkov detector, which gives an aver- 
age pion suppression factor of the order of 5 • 10 -3 . This al- 
lows to have good identification with relatively loose cuts. 
In the case of PANDA, the quartz-based Cerenkov detec- 
tor (DIRC) will provide a limited electron pion discrim- 
ination, mostly at low momentum (below 1 GeV/c): the 
necessity to have Very Tight cuts to eliminate the back- 
ground reduces the electron efficiency. 

The FFs measurement of BABAR [21] is indirect, as 
it is based on initial state radiation. The ISR correction 
factor due to hard photon emission has to be disentangled, 
and the angular dependence of the electric and magnetic 
terms is derived from elaborated simulations. The thor- 
ough study of the contributing reactions is described in 
ref. [23]. The final reconstruction efficiency is ~ 17%. 
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467SQ + 118 


+ S21 




0.008 


21 ± 01 


8.2 





2832 ± 28 


1127 ±82 


-45 ± 72 


0.035 


0.40 ±0.03 




2% 


2859 ± 28 


1057 ± 82 


45 ± 72 


0.035 


0.37 ±0.03 




5% 


2857 ± 28 


1065 ± 82 


163 ± 72 


0.036 


0.37 ±0.03 




20% 


2863 ± 28 


1042 ± 82 


769 ± 72 


0.035 


0.36 ±0.03 


13.84 





85 ±5 


39 ±19 


5± 14 


0.385 


0.47 ±0.23 




2% 


86 ±5 


41 ± 19 


12 ± 14 


0.385 


0.47 ±0.23 




5% 


86 ±5 


41 ± 19 


21 ± 14 


0.385 


0.47 ±0.23 




20% 


84 ±5 


49 ± 19 


65 ± 14 


0.382 


0.57 ±0.23 



Table 3. Results from the fit of the angular distributions, for different TPE contributions, according to eq. (|17[) . for q 2 = 5.4, 
8.2, and 13.84 (GeV/c) 2 . The first line at each q 2 corresponds to the one photon approximation. 



Accelerator 


FAIR 


CERN-LEAR 


SLAC-PEP II 


FERMILAB 


BEPC 


Experiment 


PANDA (Sim) 


PS 170 


BABAR 


E835 


BES II 


Reaction 


p + p — > e~ + e + 7 


p + p —¥ e~ + e + 7 


e~±e + — y p + p + -y 


p + p — > e~ + e + 7 


e~ + e + — > p + p 


q 2 [GeV/c] 2 


5 - 28 


3.52 - 4.18 


3.5 - 20 


8.84 - 18.4 


4 - 9.4 


£ [cm~ 2 s- L ] 


2 ■ 10^ 


3 ■ 10 3U 


3 ■ 10'" 


2 • 10' 11 


< lO^ 11 


^Beam 


10 11 p/s 


3 • 10 b p/s 




5 • 10 11 p/s 




Target 


pellets or jet 


LH2 


collider 


gas jet 


collider 


| cos 6>] 


<0.8 


< 0.8 


< 1 


< 0.62 


<0.8 


Efficiency 


40% - 10% 


~ 10 % 


17% 


67% 


~ 50% 


B/S 


< 1% 


< 5% 


< 5% 


< 2% 


1.5% - 7.8% 



Table 4. Compared characteristics of TL FFs experiments. 



7 Conclusion 

Feasibility studies of measuring proton TL FFs at PANDA 
(FAIR) have been presented. Realistic Monte Carlo simu- 
lations, which take into account the geometry, the mate- 
rial budget and the performance of the future detector, as 
well as tracking efficiency and particle identification have 
been performed. Background reactions have been studied, 
with particular attention to two body hadron production. 
The results show that, applying combined PID criteria 
and kinematical constraints, it is possible to reach a back- 
ground/signal ratio of the order of 0.1%, which is sufficient 
to ensure a clean identification of e~e + pairs correspond- 
ing to the reaction of interest ((T|) . 

We have also shown that the reaction p + p —> e~ ± e + 
at PANDA will be sensitive to a contribution of TPE of 
the order of 5% or more with statistical significance of 
about two sigma. Note also that systematical errors have 
not been taken into account. This study profits of one 
of the main advantages of FFs measurements in the TL 
region: the angular distribution of the produced electron 
in one setting contains all the useful information, allowing 
to extract the true form factors and the TPE contribution 
as well. 

The cumulated statistics, assuming four months data 
taking at the nominal luminosity for each q 2 value, will 
give precise information on the proton electric and mag- 
netic FFs, in a wide q 2 range. The precision of the ratio of 



the moduli of the electric and magnetic form factors will 
be of the order of few percent, in the overlapping region 
with the data from BABAR, which display errors of the 
order of 40%. The ratio of the electric to magnetic FF 
will be measurable until q 2 ~ 14 (GeV/c) 2 , with an error 
comparable to the existing data taken at much lower q 2 . 
Above this value it is still possible to extract a general- 
ized form factor from the total cross section and test its 
asymptotic and analytic properties. 

The measurement of the cross section allows to access 
the FFs moduli. In order to determine independently the 
real and imaginary parts, as well the relative phase, po- 
larization observables are necessary. The possibilities of 
having a polarized antiproton beam and/or a polarized 
proton target are under study. 
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